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Abstract 
Wholesale  provisioning  models  are  considered.  The  performance  measures 
used  are:  supply  material  availability,  mean  supply  response  time,  and 
availability.  These  measures  are  optimized  subject  to  a  budget  constraint. 
The  performance  constrained  budget  minimization  problems  are  also  solved  for 
the  optimal  allocation.  The  basic  method  used  is  dynamic  programming. 


I.      Introduction 
A.      Overview 

The  purpose  of  this    report   is   to  discuss   solution  methods   for  several 
problems   arising  in   inventory  provisioning.      Two  related  types   of  problems   are 
considered: 

1)  optimization  of  a  performance  measure  subject  to  a   budget 
constrai  nt 

2)  minimization  of   cost  subject  to  a   constraint   on  performance. 
The  performance  measures   considered  are  Supply  Material   Availability    (SMA), 
Mean  Supply  Response  Time    (MSRT),    and   "Pseudo-Availability"    (PA).      The  basic 
approach   to  these  problems   is   dynamic  programming. 

First   the  problems   are  formulated  and  solved  using  one   recursive 
technique.      Then  a  more  efficient   recursion   is  presented  and  discussed   in  the 
context   of  maximizing  MSRT   subject   to  a   budget   constraint. 

These  same  budget   constrained  problems   are  formulated,   discussed  and 
solved  using  a  marginal   analysis   approach    in   [2]  and   [3].     That  method, 
although  fast,   does   not   guarantee  that  optimal   solutions   are  obtained.     This 
report  provides   a  method  of  obtaining  optimal   solutions   and  thus  provides   a 
means   of  evaluating  heuristic  methods.      In  addition,   this   report  shows   how   the 
performance  constrained  budget  minimization  problems   can  be  solved  directly. 
By  contrast,   reference   [2]  addresses   this  problem  using  generalized  Lagrange 
Multipliers   or   by   solving  the  budget   constrained  performance  problem 
repeatedly   for  various   budget   levels. 

The  second  recursion  presented  here  is   considerably   faster  than  the  first 
and  it   also  guarantees   an  optimal   solution   in  some   cases.      It   is   competitive 
in  speed  with   the  marginal   analysis   method  for  small   and  medium-sized 
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problems,    but   it   is    inefficient   for  problems  with   a   large  number  of   items   or 
large  budget   values.      Its   virtue  lies   in  its  ability  to  get  exact  solutions  to 
medium-sized  problems. 

B.     Problem  Formulations ,   budget   constrained 

This   report  considers  three  specific  budget   constrained  optimization 
problems   arising  in   inventory  provisioning.      These  problems   are  formulated  and 
discussed   in  detail   elsewhere   [2]  and  will    be  stated  here  without  extensive 
explanation. 

Foil  owing   [2]  we   let 

n  =  the  total   number  of   items   considered  for  provisioning 
C,    =  the  unit   cost  of  item     i 
E.    =  the  essentiality  code  for  item  type     i 
A.    =  the  demand  rate  for   item     i 
T.    =  the  procurement   leadtime  for  item     i 
S.    =  the  number  of   items   of  type     i     provided   (the  decision 
variables ) 
Z,(S.)    =  the  performance  measure  for  item     i     when     S.      units   are 
stocked 

D.(S.)    =  Z.(S.x1v    -  Z.(S.) 

ii  i      l +1)  l  v    i  ' 

p.(x.)  =  probability  that  demand  for  item  i   is  x.  during  the 

provisioning  interval 
P.j  (x.j )  =  cumulative  probability  of  x-  or  fewer  demands  during 

a  provisioning  interval 
MTTR.  =  mean  time  to  repair  or  replace  item  i 
MTBF .  =  mean  time  between  failures  =  l/X. 
MSRT . (S . )  =  mean  supply  response  time  when  S.  units  are  stocked. 


The  three  budget   constrained  problems   considered  are: 

al)     maximize  Supply  Material   Availability    (SMA),   defined  as 


SMA(S 


1) 


,,...,Sn)    =     I     E.    X,    T.    Zr'(S.)/   T     E.    A.    T. 
1  n'        .t->      i      i      i      i  i      -S      l      i      i 

i=l  i  =1 


where 


Z 


\^(S.)    =   (1   -  p.(S.))   +   (Si    -   Xi    T.)(l   -   Pi(Si))/Xi    T. 


bl)     minimize  Mean  Supply  Response  Time   (MSRT),   defined  as 


MSRT(S,,...,S    )    =      )      E.    X.    T.    1)C'{S.)/   5"     E.    A.    T. 
v    1'  n;        . ti      i      i      i      i  i      •    i      i      i      i 

i=l  i=l 


where 


(2)  Si(Si+1 

ZUJ(S.)    =   (1   -  P.   S.        T./2  -  S./X.    +  -W — 
1  1  1       1  1  11  ?  \       T 

i      i 


+  p.(S.)(X.    T.    -   S.)/2X.    . 

Ki      i        i      i  i  i 


cl)     maximize  Pseudo-Availability    (PA),   defined  as 


rA (b ,  ,  .  .  .  »->.-  J 


z<3)(s.) 
1-1    1       ' 


where 


z<3)(s  ) 


MTBF./(MTBF.    +  MTTR .    +  Z:    y(S. ) 
i  i  i  ii 


w(2) 


In  each  of  the  budget  constrained  problems   above  the  objective  is   to 
be  optimized  by   selection  of     S.,...,S       subject  to  the  constraints 


I     Ci    S      <  B    , 
1=1     1      ] 

where  B  is  the  specified  budget  level. 


S ,    >_  U     i  nteger 


C.     Problem  Formulations,   Performance  Constrained 

In  addition  to  the  three  problems   just  stated,  we  consider  three 

related  problems   in  which  the  cost   is   to  be  minimized  subject  to  a   constraint 

on  performance. 

n 
a2)     mi  n  I     C.   S 


1=1 


i      i 


s.t.  SMA(S,  ,...  ,S    )    >  SMA 

1'        '   n     — 


S  .   >_  0     i  nteger. 


b2)     min  >P     C.   S. 

i=l     ]      1 

s.t.         MSRT(S1,...  ,S    )    <_  MSRT 

S  .    >_  0     integer. 


c2)     min  Y     C .   S. 

1-1     1      ] 

s.t.  PA(S1>... ,SR)    >  PA 


S  .    >_  0     integer, 


II.  Solution  Method  and  Examples 

A.  Dynamic  Programming  Approach 

The  computer  program  used  to  solve  these  problems  is  DP4,  a  yeneral 
purpose  program  for  performing  dynamic  programming  tabular  computations.  Here 
we  will  describe  the  general  nature  of  that  program  and  those  elements 
required  to  tailor  it  for  use  in  the  problems  considered  in  this  report. 

The  DP4  program  deals  with  a  problem  consisting  of  n  related  stages 
each  of  which  is  characterized  as  shown  in  figure  1. 

d. 


X 

1 


I 


T 


x. 


i-l 


t. (x. ,d.  ) 


r.(X.  ,d.  ) 
Figure   1.     A   single-stage  decision  problem. 


In  figure   1 


X-  is  the  "state"  variable 

d.  is  the  decision  variable 

r.  is  the  stage  return  function 

t,  is  the  stage  transformation  function. 

In  the  overall  problem  consisting  of  n  stages  the  output  state  variable  from 
stage  i  ,  namely  X.  ,  ,  is  the  input  to  stage  i-l  .  Thus,  the  n  stage 
problem  can  be  pictured  as  in  figure  2. 


'n-1 


X. 


d. 


'1 


Figure  2.      n-stage  decision  problem. 
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The  state  variable  X.  can  easily  be  understood  in  the  context  of 
the  budget  constrained  problems  as  the  amount  of  resource  (money)  remaining  to 
be  allocated  to  stages  i,  i-l,...,l  .  The  stages,  of  course,  correspond  to 
the  items  in  the  inventory  problems. 

At  each  stage  i  a  decision  d,  must  be  made.  The  decision  has  two 

effects.  First,  it  yields  a  return  r.  ,  the  performance  measure  associated 

with  the  current  item.  Second,  it  yields  a  value  of  X._,  which  serves  as 

the  input  to  the  remainder  of  the  decision  process.  The  decision  d.  must  be 

made  with  consideration  both  for  the  immediate  return  r.  and  the  future 

i 

state  X._,  .  The  overall  problem  is  to  make  the  series  of  decisions 
d  ,...,d,  to  optimize  some  function  of  the  individual  stage  returns. 

In  the  problem  (al),  where  S.  is  the  decision  variable,  we  can  let 
the  return  functions  be 


r.(X.,  S.)  =  E.  X.  T.  l[l\s.)/   I   E.  X.  T.    i  =  l,...,n 

i  =1 

and  the  stage  transformation  functions  be 


X  =  B 

n 

Xi-1  =  W  V  =  Xi  "  Ci  Si      '   =   1"">n  ' 

The  overall    return  function   is   the  sum  of  the  individual   return 
functions.     Namely, 


n 
SMA(S,,...,S    )    =     I     r.(X.,   S.)    . 


The  object   remains  to  select     S,,...,S       to  optimize  this   return. 


We  let  fi(xi)  =  the  optimal  total  return  from  stages  i,  i-l,...,l 

given  that  we  enter  stage  i  with  state  variable  X. 

i 

Then  we  can  write  the   recursive  equations   for  this   optimization  as 

i 

s.t.  X.        =   X.    -  C.    S. 

i-l  i  l      l 

and  0   <   S.    <   X./C. 

-     l   —     ii 

and  S.    =   integer 

for     i    =  2 n  . 

The  equation   at  stage   1    is 

f  1(X1 )    =  max   ^(Xp   S^ 
Sl 

s.t.  0   <_  S,    <_  X1/C1 


an 


d  S,    integer 


The  program  DP4  performs   this   optimization  provided  the  user  supplies 
the  following  subroutines   and  data. 
Required  subroutines 

1.  STGRET   -  this   subroutine  defines   the  function     r.(X.,   d.) 

ii   i 

2.  TRANFM  -  defines  the  stage  transformation  function  t.(X.,  d.) 

tit 

3.  DLIMIT  -  defines   the   range  of  decision   values     d.     which   can  be 

considered  for  the  particular  value  of     X.      under 
consideration 

4.  STORE     -  allows  the  input  of   constants   to  be  used  in  the  other 

subroutines. 


Required  Data 

1.  n  -  the  number  of  stages 

2.  For  each  stage  i 

XLOW   -  the  lowest  value  of  X.  to  consider 
XHIGH  -  the  highest  value  of  X.  to  consider 
DELX   -  the  increment  for  X. 
XMODE  -  tells  whether  to  maximize  or  minimize 

XSTAGE  -  tells  how  this  stage  return  relates  to  lower  numbered  stage 
returns  (sum,  product). 

The  methodology  is  essentially  the  same  for  the  performance  con- 
strained problems.  There  the  return  functions  r-j(X-,  d.  )  =  c.d.  .  The 
state  variable  X.  is  interpreted  as  the  portion  of  the  performance  measure 
to  be  attributed  to  stages  l,...,i  .  The  stage  transformation  functions  in 
problem  (a2)  and  (b2)  are 

xi-i  -  xi  -  W  • 

In  problem  (c2)  the  stage  transformation  is 

Vi =  yz(3)(s.) . 

The  program  DP4  and  the  subroutines  are  shown  in  Appendix  A.  The 
subroutines  are  written  to  solve  any  of  the  problems  (al),  (bl),  (cl)  or  (a2), 
(b2),  (c2).  Thus  they  involve  complications  not  needed  for  solving  just  one 
of  these  problems. 
B.  Examples 

Several  example  problems  were  solved  to  illustrate  the  approach  discussed 
here.  All  of  the  problems  involved  n  =  10  items  and  all  used  the  data  shown 
in  table  1 . 


Data 


n 

A 

Time 

Cost 

MTTR 

E 

1 

5.0 

1.0 

1.0 

.0137 

1.0 

2 

2.0 

1.0 

2.0 

.0274 

1.0 

3 

3.0 

1.0 

b.O 

.0137 

1.0 

4 

b.O 

1.0 

10.0 

.0822 

1.0 

5 

10.0 

1.0 

20.0 

.0274 

1.0 

6 

2b. 0 

1.0 

b.O 

.0027 

1.0 

7 

1.0 

1.0 

1.0 

.00b4 

1.0 

8 

1.0 

1.0 

100,0 

.0411 

3.0 

9 

O.b 

1.0 

bO.O 

.0082 

1.0 

10 

2.0 

1.0 

10.0 

.1370 

3.0 

Table  1:  Data  for  Examples 


1.  The  budget  constrained  problems. 

Tables  2,  3,  and  4  summarize  the  solutions  for  the  example  problems  (al), 
(bl),  and  (cl).  These  are  the  budget  constrained  problems. 


B  =  Budget  = 

300 

295 

290 

285 

280 

max  SMA  = 

.750654 

.746466 

.741160 

.736563 

.732964 

decision  S-.  = 

4 

4 

4 

4 

4 

S2  = 

4 

4 

4 

4 

4 

S3  = 

4 

4 

4 

4 

4 

S4  = 

5 

5 

5 

4 

5 

S5  = 

2 

2 

2 

2 

1 

S6  = 

29 

23 

27 

28 

29 

S7  " 

3 

3 

3 

3 

3 

S8  = 

0 

0 

0 

0 

0 

S9  = 

0 

0 

0 

0 

0 

S10  = 

3 

3 

3 

3 

3 

Table  2:  Solutions  to  example  problem  (al) 
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B 

= 

3UU 

295 

29U 

286 

280 

min  MSRT 

= 

.0838231 

.085962b 

.0878466 

.0900560 

.0927486 

decision 

Sl 

2 

2 

2 

2 

2 

S2 

3 

3 

3 

3 

3 

S3 

3 

4 

3 

3 

3 

S4 

4 

4 

4 

4 

4 

S5 

5 

4 

4 

4 

4 

S6 

21 

23 

23 

22 

21 

S7 

2 

2 

2 

2 

2 

S8 

0 

0 

0 

0 

0 

S9 

0 

0 

0 

0 

0 

S10 

3 

3 

3 

3 

3 

Table   3.      Solutions   to  example  problem   (bl) 


]] 


B  = 

300 

29b 

290 

285 

280 

max  PA  = 

.0403727 

.0387421 

.03700b 

.0365446 

.0341911 

decision  S, 

3 

3 

3 

3 

3 

S2 

4 

4 

4 

4 

4 

S3 

5 

4 

5 

b 

b 

s4 

b 

b 

b 

b 

b 

Sb 

3 

3 

2 

2 

2 

S6 

26 

26 

26 

27 

26 

S7 

4 

4 

4 

4 

4 

S8 

0 

0 

0 

0 

0 

S9 

0 

0 

0 

0 

0 

S10 

2 

2 

3 

2 

2 

Table  4.  Solutions  to  example  problem  (cl) 

2.  The  performance  constrained  problems. 

The  related  performance  constrained  problems  were  also  solved  for  illus- 
tration. For  example,  the  problem  (a2)  was  solved  using  the  data  from  Table  1 
with  the  restriction  that  SMA  >_  .732964  .  The  solution  to  that  problem  is  the 
same  as  the  solution  shown  in  the  last  column  of  Table  2  since  the  value 
.732964  is  the  (largest)  value  of  SMA  corresponding  to  a  budget  of  280. 
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III.     Modifications   to  Basic  Method 
A.     An  Alternative  Recursion 

An  alternative  and  more  efficient   approach   is   available  for  the  budget 
constrained  problems. 

The  approach   is   based  on  a  different   recursion  from  that  used  in 

Section   2  and   is    very   similar  to  an  approach  used  to  solve  the   "cargo   loading 

problem".      See  for  example  Dreyfus   [1]. 

The  cargo   loading  problem,    stated  as   a  maximization,    is: 

N 
max         J     v  .d  . 

j-i  J  J 

s.t.        ZCjdj<B 

d  ■   >_  0   integer  . 

Although  many  methods   are  available   for  solving  this   problem,   the  one  of 
interest  to  us   is   based  on   the  following  recursion 


f(b)   =  max         {v.   +  f(b-c)} 

je{l,...,N}      J  J 


where  f(b)  is  the  optimal  total  return  that  can  be  obtained  when  a  budget 

of  b  is  available. 

To  illustrate  this  method  consider  the  data  in  Table  5. 


v . 


1 


3  =  N 


=  1U 


Table  b.  Data  for  example  using  alternative  recursion 


The  computation  proceeds  with  increasing  values  of  b  until  b  =  B  is 
reached.  The  process  can  be  viewed  as  shown  in  figure  3  where  a  template 
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representing  the  available  items  is  placed  over  the  budget  value  of  current 
interest.  The  template  points  back  to  previously  determined  optimal  solution. 
Each  of  these  previous  solution  is  considered  for  updating  by  including  one 
more  item  of  the  type  inducted  by  the  template.  The  updated  solutions  are 
compared  and  the  best  selected  as  the  solution  for  the  current  value  of  b  . 
The  illustration  shows  the  template  at  the  budget  value  of  7.  The  optimal 
solutions  for  b  =  0,1,. ..,6  have  already  been  computed.  The  comparison  at 
B  =  7  is  among  the  solution  at  6  with  an  additional  item  1  for  a  total  return 
of  9,  the  solution  at  4  with  an  additional  item  of  type  2  for  a  total  return 
of  10,  and  the  solution  at  3  with  an  additional  item  of  type  3  for  a  return 
of  10.  Either  of  the  last  two  is  chosen  and  recorded  as  an   optimal  solution 
at  7. 


d,(b) 

di(b) 
do(b) 


f(B)  =  0 


10 


10 


*an  alternative  optimal 
solution  at  b  =  6  is 
d,  =  0 


1 


d2  =  2 


d3  =  0 


Figure  3:  Illustration  of  solution  method  using  the  alternative  recursion. 

A  very  simple  modification  of  this  procedure  can  be  used  to  solve  the 
problems  of  the  type  discussed  in  this  report.  We  consider 


max  I     r.(d.) 
j=l  J  J 


s.t.   7  c,  •  d.  < 

.I-'.  1  1  — 


J-l 


J   J 


d  .  >_  0  integer 


l-i 


where  the  return  functions  ri(d-)  are  concave.  The  procedure  below  is  not 

guaranteed  to  give  the  optimal  solution  for  all  r.(d.)  but  is  guaranteed 

if  the  r-;(d-j)  are  points  on  a  concave  function. 

In  this  case  we  will  represent  the  return  functions  r.(d.)  as  the 

J     J 

sum  of  the  marginal    values   of  additional   items   of  type     j   . 


r.(d)   =  +     I     m,(1) 


i=l 


Thus  r.(2)  =  m.(U)  +  m.(l)  +  m.(2)  .  These  marginal  values  form  a 
sequence  with  the  properties  that 


and 


m.(i  )  >  0 
J 


m.  (i  )  >  m.  (k  )     i  <  k 
J      J 


The  same  algorithm  as  before  was  applied  with  the  modification  that  the  value 
term  v-  ,  which  was  formerly  constant,  is  replaced  by  m-(i)  for  the 
appropriate  value  of  i  .  The  program  which  implements  this  algorithm  was 
called  RECUR. 

It  should  be  noted  that  the  discussion  above  treats  the  constraint  as 
an  inequality,  but  the  function  f(b)  in  this  section  is  computed  for  the 
constraint 

N 

T  c  .x .  =  b  . 

J=l  J  J 

For  this   reason  we  may  have   in  a  maximization  problem 

f(bx)    >  f(b2) 


lb 


although  b,  <  b~  .  That  is,  if  we  require  the  equality  to  be  met  exactly, 
it  is  not  necessarily  true  that  a  larger  budget  is  better.  The  program  prints 
the  values  of  f(b)  for  several  values  of  b  so  the  optimal  value  of  b  <_  B 
can  be  found  visually. 
B.  Minimum  Orders 

1)  RECMOD 

A  modification  was  made  to  the  program  RECUR  to  permit  the  user 
to  specify  minimum  packaging  quantities  of  each  item.  That  is,  item  i   is 
assumed  to  be  packaged  with  q-  items  per  package.  The  provisioning  can 
select  only  whole  packages  of  each  item.  This  modification  resulted  in  the 
program  RECMOD  given  in  the  Appendix. 

2)  Example 

To  illustrate  the  RECMOD  program,  consider  the  example  problem 
(a2)  solved  previously.  The  optimal  solution  for  budget  B  =  300  is  repeated 
in  Table  2  for  the  case  in  which  all  q.  =  1  . 


n 

1 

3 

3 

4 

5 

b 

7 

8 

9 

10 

Objective 

^ 

1 

1 

i 

1 

1 

1 

1 

1 

1 

1 

di 

2 

3 

3 

4 

5 

21 

3 

U 

0 

3 

.0838231 

«1 

2 

1 

1 

5 

i 

1 

1 

3 

1 

1 

di 

2 

3 

3 

5 

4 

33 

3 

U 

0 

3 

.0847410 

Table  b.  Solution  to  Example  using  RECMOD. 

The  optimal  solution  is  also  shown  for  a  modified  problem  in  which  not  all  q. 
are  equal  to  1  . 
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C.  Discussion 

The  modified  recursion  just  discussed  has  been  implemented  for  the  prob- 
lem (a2)  which  is  to  minimize  MSRT  subject  to  a  budget  constraint.  The  method 
can  also  be  applied  to  the  other  budget  constrained  problems,  but  this  has  not 
yet  been  done.  All  that  is  required  is  to  modify  the  program  to  compute  SMA 
or  PA  instead  of  MSRT  and  to  maximize  instead  of  minimize. 

There  is  a  difficulty  in  extending  this  method  to  problems  in  which  the 
item  costs  are  arbitrary  values.  The  method  is  very  effective  when  the  costs 
are  all  integer  and  can  be  scaled  so  that  the  smallest  cost  is  1  .  If 
arbitrary  costs  are   allowed,  the  algorithm  can  become  ineffective  for  all 
except  small  values  of  B  .  Consider  for  example  the  costs  of  $1.00,  $1.21, 
$1.27  for  three  items.  Let  the  budget  be  52b. UO.  The  problem  could  be  solved 
by  scaling  the  costs  to  be  100,  120,  and  12b  and  the  budget  to  be  2b00,  but 
then  too  many  values  of  b  must  be  considered  when  many  of  them  are  not  pos- 
sible. Alternatively  the  algorithm  can  step  to  the  "next  possible  value"  and 
will  consider  the  folowing  sequence  of  values 

b  =  100,  121,  127,  200,  221,  227,  242,  248,  2b4,  300,  321,  327,...  . 
As  the  process  continues,  depending  on  the  relative  values  of  the  costs,  the 
sequence  becomes  more  dense  and  may  eventually  include  all  possible  values  of 
b  .  This  is  ineffective  and  cumbersome  for  large  values  of  B  . 

It  is  also  not  possible  to  apply  the  modified  recursion  for  the  perform- 
ance constrained  problems.  On  the  other  hand,  the  budget  constrained  problems 
are  solved  very   rapidly  and  the  relationship  between  cost  and  performance  can 
easily  be  determined  from  solving  the  budget  constrained  problem  using  the 
first  recursion.   In  fact,  if  MSRT  is  minimized  for  a  budget  of  B  ,  the 
solution  is  actually  obtained  for  all  values  of  b  up  to  and  including  B  . 
Those  results  reveal  the  relationship  between  performance  and  budget. 
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APPENDIX-PROGRAM  LISTINGS 


FILE:    DPFILE 


exe: 


Al       NAVAL    POSTGRADUATE    SCHOOL 


FILEDEF 

07 

DISK 

{PERM 

FILEDEF 

08 

DISK 

RUSS10     DATA     (PERM 

FILEDEF 

06 

TERM 

(HECL     133    RECFM    FB 

FILEDEF 

02 

DISK 

(PERM 

PE*M 


FILE: 

FILE 

10 

NEW 

1 

10 

300. 

10 

CLD 
2  95. 

10 

CLD 
290. 

10 

OLD 
285. 

10 

OLD 
280. 

10 

CLD 
275. 

FT07F301    Al       NAVAL    POSTGRADUATE    SCHOOL 


0. 


3  00. 

300. 

2  9  5. 
2  90. 
235. 
280. 
2  7  5. 


i. 


MIN 


SUM 


FILE:    RUSS10 


DAT4 


Al       NAVAL    POSTGRADUATE    SCHOOL 


2.  2. 

3.  5. 
5.                    13. 

10.  20. 

25.  5. 

1.  1. 

1.  100.  3 
.5                   53.  1 

2.  13.  3w 
10       303(06, CON)3C0ST    1SMA    2MSRT    3A^IL   - 

21151113 
FOR    RECUR    THE    LAST    TWO    LIME    MUST    BE    PLACED 
CURRENT    LINE.       IT     IS    FOR    3P4 


.0137 
.327^ 
.0137 
.0822 
.027^ 
.0027 
.005^ 
.0411 
.0082 
.1370 
DATA    ACEB3  =  LAM, COST,  ESS, T,M? 

FIRST    INSTEAD    OF    THE 
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C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


DATA    DESCRIPTION     FOR    CARD    NUMBER     1 


COLUMNS 

.    UUSTIFY    , 

.    VARI BLE    . 
»    NAME 

.    MEANING 

1-5 

,    C3L    5 

.    MMAX 

.    THE     NUMBER    OF    STAGES    IN    THE     PROBLEM 

6-10       , 

,    COL    10 

.MTAPE 

,     LOGICAL    TAPE    NO    OF     THE    MAIN    TAPE 

.     IF     THIS     IS    LEFT    BLANK,     THE    COMPUTER 

.     WILL    USE    STANDARD    SCRATCH   TAPE    3 

11   -     15    , 

►    COL    15 

.    SOLVE 

.    =      OLD          IF    THE    OPTIMAL    DECISION 

FUNCTIONS    HAVE     ALREADY     BEEN 
CALCULATED    AND    ARE     LOADED 
ON    LOGICAL     TAPE    MTAPE 

►    =       NEW          IF    THE    OPTIMAL     DECISION 
FUNCTIONS     MUST     BE 
CALCULATED    AND    STORED    ON 
LOGICAL    TAPE    MTAPE    BEFORE 
SOLVING    THE    PROBLEM 

DATA    DESCRI3TI3N     FOR    THE    NEXT    G^OUP   OF    CARDS 

IF    SOLVE    =    NEW       ON    CARD    NUM3ER    1,      THEN    EACH    OF    THE    NMAX    STAGES 
MUST     BE    DEC^IBED     BY    THE    CARDS    DISCUSSED    BELOW.    YOU    CAN    USE    1    CARD 
PER    STAGE    OR    YOU    CAN    MAKE     1    CARD    DESCRIBE    MAMY    ADJACENT    STAGES    IF 
THEY     ARE    SIMILAR.    THE     STAGES    MUST     BE    DESCRIBED    I N*    NUMERICAL    ORDER 
STARTING    WITH    STAGE    1     (I.E.    STAGE     1,     THEN    2,     ...,     THEN    NMAXi 
OMIT    THIS    PACK    OF    CARDS    WHICH    DESCRIBE    THE    STAGES    IF    SOLVE    =    OLD 


COLUMNS    , 
1-5 
6-10       , 

11    -    20    . 
21    -    30    . 

»    JUSTIFY    , 
.    COL    5 

■    COL    10       , 

,    ANY 
ANY 

,    VARIBLE    , 
,    NAME 

.NSTAGE       , 
,    NDITTO       ! 

,    XLOW 
,    XHIGH 

■     ME  4  N I NG 

LOWEST    NUMBERED    STAGE    FOR    WHICH 
,    HI  S    CARD    APPLIES 

,    HIGHEST    NUMBERED    STAGE    FOR    WHICH 
►    THIS    CARD    APPLIES.     IF    THESE    COLUMNS 
ARE     LEFT    BLANK,     THEN    NDITTO     WILL    BE 
,    TA<EN    AS    EQUAL    TO    XSTAGE 

,    LOWEST    VALUE    OF    XM    FOR    STAGE    NSTAGE 

,    HISHEST     VALUE    OF    XN    FOR    STAGE   NSTAGE 

31    -    40    , 

,     ANY 

,    DELX 

,     INCREMENT    IN   XN    FOR    STAGE   NSTAGE 

41    -     46    , 

,    COL    46       . 

►    XMODE 

,    =         MIN       IF    STAGE    NSTAGE    IS    TO    BE 
MINIMIZED 

,    =         MAX       IF    STAGE    NSTAGE    IS    TO    BE 
MAXIMIZED 

47   -     52    , 

►    COL    52 

,    XSTAGE       , 

,    =         SUM       IF    THE    COMPOSITION 

OPERATOR    BETWEEN    STAGES 
NSTAGE    AND    NSTAGE-1    IS 
ADDITION 

.    =       MULT       IF    THE    COMPOSITION 

OPERATOR    BETWEEN   STAGES 
NSTAGE    AND    NSTAGE-1    IS 
MULTIPLICATION 

.    =41 NMAX        IF    THE    COMPOSITION 

OPERATOR    BETWEEN    STAGES 
NSTAGE    AND    NSTAGE-1    IS    TO 
MINIMIZE    THE    MAXIMUM 
INDIVIDUAL    STAGE    RETURN 
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C 

c 
c 
c 

c 
c 
c 
c 

C 

C 
C 
C 
C 

c 

C 

c 

C 

c 

c 
c 
c 
c 
c 
c 

C 

c 
c 
c 

C 

c 

c 
c 

C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

C 

c 
c 
c 
c 
c 
c 
c 
c 
c 

C 
C 

c 
c 

c 
c 

c 
c 
c 


►    =M*XMIN       IF    THE    COMPOSITION 

OPERATOR    BETWEEN   STAGES 
NSTAGE    AMD    NSTAGE-1    IS    TO 
MAXIMIZE    THE    MINIMUM 
INDIVIDUAL    STAGE    RETURN 

,    =                      THIS    VARIABLE    IS    NOT   USED 
WHEN    NSTAGE=1,     THEREFORE 
COLUMNS    *7    TO    52    ARE 

IGNORED    WHEN    NSTAGE=1    AND 
MAY    BE    LEF  T    BLANK 
HOWEVER    IF    NDITTO    IS 
GREATER   THAN    It    THEN   XSTAGE 
MUST    BE    SPECIFIED    IN 
ACCORDANCE    WITH    THE 
COMPOSITOR    OPERATOR   FOR    THE 
REST    OF    STAGES    DESCRIBED   ON 
THIS    CARD 

53    -    55    , 

.    COL    55       , 

.MODES 

,    =      1               IF    THE    FIBONACCI    SEARCH    IS 
TO    BE    USED 

,    =      2               IF    THE    FIBONACCI    SEARCH    IS 
NOT    TO    BE    USED 

DATA    DESCRIPTION    OF    LAST    CARD    F3*     PROBLEM 

THE    OPTIMUM    DECISIONS    ARE    TO    BE   PRINTED    OUT    FOR    A    PROCESS    WITH 
NMAX    STAGES,     FURTHER,     FN(XM     IS   TO    BE     OPTIMIZED    WITH    RESPECT    TO  XN 
WHEN    THE    PROGRAM    COMPLETES    THE    REQUESTED    PRINTOUT    IT    WILL     RETURN 
TO    THE    BEGINNING    AND    START    A    NEW     PROBLEM    IF    THERE     IS    DATA,    IF    NO 
DATA,     IT    WILL    EXI T 


COLUMNS    , 
1    -    15      , 
16    -    30    , 

.    JUSTIFY    , 
,    ANY 
,    ANY 

,    VARIBLE    . 
,    NAME 

,    XNl 
,    XN2 

,     ME  A  N  I NG 

,    AT    STAGE    NMAX,     OPTIMIZE    FN(XN)     FCR 
.    XSJ    NOT    LESS    THAN    XNl 

,     AT     STAGE    NMAX,     OPTIMIZE    FN(XN)     FCR 
,    XN    NOT    GREATER    THAR    XN2 

DESCRIPTION    OF    THE    SUBROUTINE    T^ANFM 

THE    MAIN    PROGRAM    TRANSMITS    XN,DN,     AND     NUMBER     (THE    SEQUENCE    NUMBER 
OF    THE    PROBLEM*     TO   THE    SUBROUTINE     AND    THE    SUBROUTINE    CALCULATES 
YN    (THE    OUTPJT    OF    THE     STAGE    I.E.     YN=XN-1) 


DESCRIPTION    3F    THE    SUBROUTINE    STGRET 

THE    MAIN    PROGRAM    TRANSMITS    XN,DN,YN,    AND    NUMBER    TO    THE 

AND    THE    SUBR3UTINE    CALCULATES    RN     (THE    STAGE    RETJRN) 


SUBROUTINE 


DESCRIPTION    3F    SUBROUTINE    STORE 

THE    MAIN    PROGRAM    CALLS    THE    SUBR3JTINE    RIGHT    AFTER    READING     THE 
FIRST    DATA    CARD    AND    DEFINING    NUM3ER.     THE    PURPOSE    OF    THIS 
SUBROUTINE    IS    TO    ALLOW    THE    STORING    OF    CONSTANTS    IN    COMMON    STORAGE 
FOR    THE    POSSIBLE    USE    OF    SUBROUTINES    STGRET    AND   TRANFM 


DATA  WHICH  WILL  BE  READ  FRCM 
THE  SUBROUTINE  STORE  MUST  BE 
AND    THE    GROUP    OF    CARDS    WHICH 


DATA    CARDS    AT    DBJECT    TIME    BY 
INSERTED    BETWEEN    THE    FIRST   DATA    CARD 
DESCRIBE    THE    INDIVIDUAL    STAGES 
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5  DESCRIPTION    DF    SU6R0JTINE    DLIMIT 

C  THE    MAIN    PRD3RAM    TRANSMITS    N.XN,     AND    NUMBER    TO   THE    SU6R0UT IME , 

C  AND    THE    SUBR3UTINE    CALCULATES    DLDW,    3HIGH,    AND    DELD 

C  WHERE    DLOW    =    LOWEST    VALUE    OF    DN    FOR    THAT    PARTICJLAR    VALUE     OF    N    AND 

C  XN 

C  DUGH    =    HIGHEST    VALUE    OF    DM        FOR     THAT    PARTICULAR    VALUE    OF    N 

C  AND    XN 

C  DELD    =    INCREMENT    IN    DN    FOR    THAT    PARTICULAR    VALUE    OF     N   ANC   XN 

C  IF    ANY    OF    THESE    VARIBLES     DO    NOT    ;HAN3E,    THEY    CAM     BE    SET    IN 

C  SUBROUTINE    STORE    AND    NOT    MENTIONED    IM    SUBROUTINE    DLIMIT 

C 

C 

C 

C  IF    ANY    SUBROJTINE    SETS    MUFF   TO   -1       THEN    THE    PR03RAM    IGNORES   THE 

C  VALUE    OF    DN    BEING    PROCESSED    AT    THAT    FOMENT,    AND    IGNORES    ALL    LARGER 

C  VALJES    OF    DN    FDR     THAT    PARTICULAR     VALUE    OF    XN    AT    THAT    PARTICULAR 

C  STAGE — THIS    FEATURE    IS    DISABLED    dHEN    USING    THE    FIBONACCI     SEARCH 

C 

C 

c 

C  THIS    PROGRAM    USES    TAPES    3    AND    4   FOR   SCRATCH 

C  NSS2    ZERO    FOR    ENTIRE    PROBLEM    SO-</ING    AND    OUTPUT 

C  NSS2     NON-ZERO    FOR    OUTPUT     FROM    PREVIOJS    PROBLEM 

DIMENSION    FN (  10001)  ,  FNMK  10301 ), 3 NQFX N( 1 0001 ) 

REAL    *    8    H0L1IN/1  HIN*  /,     H0L0L3/*       OLD'/ 

REAL    *    8    XM03E,     XSTA3E,    SOLVE 

REAL     *    8    TAB(4) 

COMMON    /CON/TAB 

COMMON    XN,N,3N, YN  ,RN, SLOW  ,XH I GH, 3 ELX,  DLOW, DH I GH, DELD, NUMBER 

COMMON    YLOW,YHIGH,DELY,NUFF,FN,JrOP,MNF«_AG,NMAX,XLAM 

COMMON    A,B,C,D,E, SUME LT. KQDElt K30E2 

EQUIVALENCE    {  FN  (  1  ) ,  FNM1(  1  )  ,DNOFXM(  11) 

WRITE(6,1000)     HOLMIN,     HOLOLD,     (TABU),     1=1,4) 
1000    FORMATdX,    646) 

MAX= 10001 

INTAPE=7 

N0UTPE=6 

CALL     S  EARCHO.  0  t  0.0 1  -  1. 0  t  0.0  t  DN3EST.BE  ST) 

NPAGE=0 

NUMBER=0 
240    READ( INTAPE,103 ,END=2000 )     NMAX,^TAPE,     SOLVE 
100    FOPMAT(2I5,A5) 

MTAP    =    MTAPE 

NDITTO=0 

NCALC=4 

IF<MTAPE-4)    405,501,405 
501    NCALC=3 

GO   TO    172 
405    IF(MTAPE)     322,322,172 
322    MTAPE=3 
172    NTAP  E=MTAPE 

REWIND    MTAPE 

NUMBER=NUMBER+1 
C  CALL    STORE    WAS    MOVED    FROM    HERE 

IF(SOLVE-HOLOLD)     301,302,301 

302  NSS2=1 
NPAGE=NPAGE+1 
WRITE(N0UTPE,152)     NPAGE 

152    FORM AT(1H1103X4HPAGEI4) 
WRITE(NOUTPE,303)     MTAPE 

303  FORMAT    (58H    THE    TABLES    OF    FN(XN)     AND    DN(XN)    WHICH    ARE    ON    LOGICAL    T 
1APEI5,35H    WILL    BE    USED   TO    SOLVE   THIS    PROBLEM) 

READ(MTAPE)    MMAX 

REWIND    MTAPE 

GO    TO    304 
301    NSS2=0 

CALL    STORE 

NPAGE=NPAGE+1 

WRITEiNOUTPE, 152)     NPAGE 

WRITE(NOUTPE,305)     MTAPE 
305    F0RMATO9H    THE    COMPUTER    IS    TO    CALCULATE    TABLES    OF     FN(XN)     AND    DN(XN 
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1),    AND    STORE    THEM    ON    LOGICAL    TAPEI5) 
REWIND    NCALC 
500    WRITE(^TAPE)     NMAX 
304    WRITE(NOUTPE, 306)     NMAX 

306   F0RMATU6H   THE    PROBLEM   HASI5t7H   STAGES) 
314    IFCNSS2)     197,109,197 

109  LINES    =    6 

WRITE(N0UTPE,2006 )    NMAX,     MTAP,    S3LVE 

2006  FORMAT(///    •     DATA    CARD       • ,    217,    47) 
121    DO    123    N=1,NMAX 

NM1=N-1 

IF(NDITTO-N)     241,132.132 
241    REA0(INTAPE,124)     NSTAGE, NDITTO, XL OW,X HIGH , DEL  X, XMODE, XSTAGE, MODES 

124  F0RMAT(2I5,3E10.6,2A6 ,13 ) 

WRITE(NOUTPE,2007 )    NSTAGE,    NDITTD,    XLOW,     XHIGH,    DELX,     XMODE, 
1  XSTAGE,    MODES 

2007  FORMAT!///    '    DATA    CARD       ',    215,    2X,    3F10.6,    2Ab,     13) 
IF(LINES-41)     30,30,31 

31    NPAGE=NPAGE+1 

WRITEtNOUTPE, 152)    NPAGE 

LINES=1 
30    LINES    =    LINES    +    9 

ITOP=( XHIGH- XLOW) /DELX +1.001 

IF(MODES-l)    600,601,600 

600  MODE  S=-l 

601  IF(NSTAGE-N)     125,245,125 

125  WRITE(N0UTPE,127) 

127    F0RMAT(23H1DATA    CARD    OUT    OF    ORDER) 
GO   TO    2003 

245  IF(NSTAGE-NDITTO)     246,126,126 

246  WRITE(N0UTPE,247)    N,N  STAGE  ,NDI  H3 

247  F0RMATC///46H    THE    DESCRIPTION    WHICH    APPEARS    BELCH     FOR     STAGE, 15, 17H 

I  APPLIES    TC    STAGE, 15, 11H    THRU    STAGE,  15, 1  OH    INCLUSIVE) 
LINES=LINES+4 

XXMODE    =    XMODE    -    HOLM  IN 
WRITE( 6,1001)     XXMODE 
1001    FORMATUX,    A6  ) 

126  IF(XMQOE-HOLMIM)     230,231,230 

231  XM00E=1.0 
WRITE(N0UTPE,232 )    N 

232  F0RMAT(//6H    STAGEI5,18H    IS    A    MINIMIZATION) 
GO   TO    520 

230    XM0DE=-1.0 

WRITEtNOUTPE, 234)     N 
234    F0RMAT(//6H    STAGE  15, 18H    IS    A    MAXIMIZATION) 

520  00    129    1=1,4 
IF(XSTAGE-TA3( I) )     129,130,129 

129  CONTINUE 

521  NNFLAG=5 
GO   TO    502 

130  NNFLAG=I 

502    IF(ITGP-MAX)     316,316,112 

316    WRITE(N0UTPE,110)     XLO  W  ,XHI  GH  ,DEL< 

110  F0RMAT(63H    THE    STATE    VARIBLE    IS    BEING     TREATED    AS    OISCRETE,    GOING    F 
1R0M    XN  =  E14.5,7H    TO   XN=E14.5,13H    IN    STEPS    OF    E14.5I 

IF(MOOES)     318,318,321 
313    WRITE(NOUTPE,510) 

510  F0RMAT(41H    TDTAL     ENUMERATICN    IS    J S ED    FOR    THIS    STAGE) 
GO   TO    320 

321    WRITE(N0UTPE,511) 

511  F0RMAT(40H    FIBONACCI    SEARCH    IS    JSED   FOR    THIS    STAGE) 
320    GO   TO    (11,12, 13, 1 4, 132 ) , NNFLAG 

11  WRITE(N0UTPE,16 )     NM1,N 

16  F0RMAT(40H    THE    COMPOSITION    OPERATOR    BETWEEN    STAGESI6,4H   ANDI6,12H 
IIS    ADDITION) 

GO    TO    132 

12  WRITE(NOUTPE,  17)     NM1,  N 

17  F0RMAT(40H    THE    COMPOSITION    OPERATOR    BETWEEN    STAGESI6,4H   ANDI6,18H 

II  S   MULTIPLICATION  ) 
GO   TO    132 

13  WRITE(N0UTPE,18)     NM1 , N 

18  F0RMAT<40H    THE    COMPOSITION    OPERATOR    BETWEEN    STAGESI6,4H   ANDI6,34H 
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IIS    TO    MAXIMIZE    THE    MINIMUM    RETURN) 
GO   TO    132 
14    WRITE(N0UTPE,19)     NM1 , N 

19    FORMAT(40H    THE    COMPOSITION    OPERATOR    BETWEEN    STA3ESI6,4H    ANDI6,34H 
IIS    TO    MINIMIZE    THE    MAXIMUM    RETURN) 
132    DO    136    I  =  l,ITOP 
134    BEST=1.0E+35*XMODE 
X  =  I-1 

XN=XLOW+X*DELX 
NUFF=1 
CALL     OLIMIT 
IF(MODES)    602,602,603 
603    CALL    SEARCH(DLOW,DHIGH,OELD,XMODE,DNBEST,BEST  ) 

GO    TO    400 
602    KTOP=( DHIGH-3L3W) /DELD+1.001 
223    00    137    J=1,KT0P 
X=J-1 

DN=DLOW+X*DELD 
C  JUMP    MAY    BE    SET    TO    1     IN    TRANFM   TO     INDICATE   AN    IMFtASIBLE    XN-1 

C  JUMP=0 

CALL     TRANFM 
CALL    STGRET 
C  XM0DE=1    FOR    M I M  ,     -1    FOR    MAX 

C  I F( JUMP. EG. 1) RN=999999*XM0DE 

IF(NUFF)     400,401,401 

401  IF(N-1  )•  504,504,503 
504   QN=RN 

GO   TO    143 
503    K=(YN-YLOW)/3ELY+1.001 

I F<K  »     212,212,209 
212    WRITE(NOUTPE,210 J    N,XN,DN,YN 

210    F0RMAT(///9H    AT    STAGE  15, 9H    WITH   KN=E15,8,8H    AND    DN=E15.8,12H    XN-1 
1EQUALSE16.8,21H    AND    IS    OUT    JF    LIMITS///) 

GO   TO    2003 
209    IF(K-JTOP)    205,206,212 

206  RNM1  =  FNMKJT3P) 
GO    TO    207 

205    X=K-1 

X=YLOW+X*DELY 

RNM1  =  FNMHK)*(FNM1(K>1)-FNM1(K))  *  (  YN-X)/DELY 

207  GO    TO     (21,22,23, 24, 248), NNFLAG 

248  WRITE(N0UTPE,249)     N 

249  F0RMAT(43H1C3MP0SITI3N   OPERATOR    MOT   DEFINED    FOR    STAGE, 15) 
GO   TO    2003 

21  QN=RN+RNM1 
GO   TO    143 

22  QN=RN*RNM1 
GO   TO    143 

23  IF(RN-RNMl)    141,141,142 

141  QN=RN 

GO   TO    143 

142  QN=PNM1 
GO    TO    143 

24  IF(RN-RNMl)     145,145,146 

145  QN=RNM1 
GO   TO    143 

146  ON=RN 

143  IF(( GN-BEST)*XMODE)  144,137,137 

144  BEST=QN 
DNBEST=DN 

137    CONTINUE 

400    WRITE(NCALC)     BEST,DNBEST 

136    CONTINUE 

REWIND    NCALC 

WRITE(NTAPE)     XL OW ,XHI GH, DELX , ITD? ,  XM3DE 

DO    404    I=i,ITOP 
404    READ(NCALC)    DUM, DNOFXN( I ) 

REWIND    NCALC 

WRITE (NTAPE)     (DNGFXN(  II)  , I  I  =  1,IT0P) 

DO   402    1=1, ITOP 

402  READ (NCALC)    FN(I),DUM 
REWIND    NCALC 
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YLOW=XLOW 
YHIGH=XHIGH 
DELY=DELX 
JTOP=ITOP 
123    CONTINUE 

WRITE(NTAPE)     ( FN ( II ) , 1 1=1 , ITOP ) 

161  ENO    FILE    MTAPE 
197    REWIND    MTAPE 

READ(MTAPE)    NMAX 

N=NMAX 

J=N-1 

IF(J)    407,407,201 
201    DO   522    1=1, J 

READ(MTAPE) 
522  READ(MTAPE) 
407    READ(MTAPE)    XLOW , XHI3H,DE LX, ITOP  ,  XMOO E 

READ (MTAPE) 

READ(MTAPE)     (  FM  (  I  )  ,  1=1 , 1  TOP  ) 

BACKSPACE    MTAPE 

BACKSPACE    MTAPE 

READ  (INTAPE, 162)     XN1,XN2 

162  F0RMAT<2E15.8) 
NPAGE=NPAGE+1 
WRITE(NQUTPE,152)     NPAGE 
WRITE(NOUTPE,2008)    XN 1 ,     XN2 

2008    FORMAT<///    ■     DATA    CA^  D       •,    2F20.3     /// ) 
IF(XMODE)     323,323,324 

324  WRITE(NOUTPE,325)     N 

325  F0RMAT(29H    THE    PROBLEM    IS    TO    MINIMIZE     AI6,14H    STAGE    PROCESS) 
GO   TO    327 

323    WRITE(NOUTPE,326)    N 

326  F0RMAT(29H    THE    PROBLEM    IS    TO    MAXIMIZE     AI6,14H    STAGE    PROCESS) 

327  WRITE(N0UTPE,32S)     XNi,XN2 

328  FORMAT    (42H    XN    IS    TO    BE    CHCOSEN   OPTIMALLY    BETWEEN    XN= 1PE1 4.5, 8H   AN 
ID    XN  =  1PE14.5) 

i71    IF(XLOW-XNl)     175,175,176 

176  WRITE(N0UTPE,177 I     XLQ^,XHIGH 

177  F0RMAT(47H0THE    PROGRAM    ONLY    HAS    INFORMATION    ON    XN'  BET  WEENE16.  8  ,4H 
1ANDE16.8) 

GO    TO    240 
175    IFUN2-XHIGH)     181,181,176 
181    IFUN1-XN2)    178,173,176 

178  I  Xl=  (XNl-XL3w)/DELX-H.001 
1X2= (XN2-XL0WJ/DELX+1.001 
BEST=1.0E+35*XM0DE 

DO    182    J=IX1, 1X2 

IF((  FN(J) -BEST >*X MODE  )    183,182,13  2 
133    BEST=FN(J) 

JSAVE=J 
132    CONTINUE 

X=JSAVE— 1 

XN=XLCW+X*DELX 

WRITE(N0UTPE,184)     XN, BEST 
184    F0RMAT(12H    OPTIMAL    XN  =  IP  E  14.  5,  1SH     OPTIMAL    RETURN=  1PE1  4.  5) 

READ  (MTAPE)     (DNOFXNd  I  ) ,  I  1=1 ,  ITD?  ) 
186    DN=DNOFXN(JSAVE  ) 

CALL  TRANFM 
204    WRITE(N0UTPE,188> 

188  F0RMAT(7H0  Nl 8X2H XN1 3X2HDN1SX4HXN-1 ) 
NN=0 

203    WRITE(N0UTPE,189)     N,XN,DN,YN 

189  FORMAT( I7,1P3E20. 5) 
NN=NN+1 

N  =  N-1 

IF(N-l)    244,193,193 
244    GO    TO    240 
193    DO    406    1=1,4 
406    BACKSPACE    MTAPE 

READ (MTAPE)    XLOW , XHIGH,DE LX, ITOP  ,  XMOOE 

READ (MTAPE)     (DNOFXNi I  I ) , I  I =1 , 1 TD3 ) 

L=( YN-XLOW)/D EL X+ 1.001 

IF(L)    214,214,215 
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214  L  =  l 

219  NP1=N+1 
WRITE(NOUTPE,210)     NP1,XN,DN,YN 
NN=NN+7 

215  IF(L-ITOP)    218,216,217 

216  DN=DNOFXN(L) 
XN=YN 

GO    TO    196 

217  L=ITOP 

GO    TO    219 

218  XN=YN 
194    X=L-1 

X=XLOW+X*DELX 

DN=DNOFXN(L)MDNOFXN(  L+l  )-ONOFXM(L)  )*  {  XN-X  )  /DELX 
196    CALL    TPANFM 

IF<NN-50)    203,203,220 

220  NPAGE=NPAGE+1 
WRITE(N0UTPE,152  )     NPAGE 
GO   TO    204 

112  WRITE(NOUTPE, 113)     MAX,ITOP 

113  F0RMAT(72H    THIS     PROGRAM    LIMITS    THE    NUMBER    OF     DISCRETE     STEPS   OF    THE 
1STATE    VARIBLE    T0I6,26H    AND    THIS   PROBLEM    REQJIRESI6) 

2002  FORMATCl','     ERROR    HALT1) 

2003  WRITE(6,2002) 
STOP 

2000  WRITE(6,2001) 

2001  F0RMAT(•1•,,       END    OF    DATA    FILE') 
STOP 

END 
CSEARCH  DIS3RETE    FIBONACCI    SEAR:H    SUBROUTINE 

SUBROUTINE    SE ARCH ( AA, BB, DELYY, XXiODE,  YYBEST, 8  EST ) 

DIMENSION    FI150I 
C 

C  OPTIMIZE    WITH    RESPECT    TO    Y    BETWEEN    AA    AND    BB    IN    STEPS    OF    DELYY 

C  STORE    OPTIMUM    Y     IN    YYBEST    AND    THE     OPTIMUM    VALUE    OF    THE    OBJECTIVE 

C  FUNCTION    IN    BEST 

C  XXM0DE=-1    FOR    MAXIMIZE 

C  XXM0DE=1       FOR    MINIMIZE 

C  SUBROUTINE    MJST    BE    INITIALIZED    3*     CALLING    IT    WITH    DELYY=-1.0    AT 

C  LEAST    CNCE    BEFORE    IT     IS    USED    FG*     A    SEARCH 

DELY=DELYY 

A=AA 

XMODE=XXMODE 
IF(DELY)     100,100,101 

100  F(1)=1»0 
F(2)=2.0 

DO   1    1=3,150 

F( I)=F( I-ll+F ( 1-2  1*1. 0 

IF(FU  J-l.OE+35  )     1,2,2 
2    I  1  =  1 

RETURN 
1    CONTINUE 

11  =  150 

RETURN 

101  YHI=B 
YLO=A 

FNO= (  YHI-YL0)/DELY+1.0 
DO    5    N=1,II 
IF{FNO-F(N))    6,6,5 

5  CONTINUE 
WRITE( 6,7) 

7    F0RMAT(38H1ERR0R,    TOO    MANY    POINTS     TO    BE    SEARCHED) 

2002  F0RMAT(•1,,,     ERROR   HALT* I 

2003  WRITE* 6,2002) 
STOP 

6  IF(N-2)     42,41 ,40 
40    Y1=F(N-2)*DELY+YL0 

CALL     FUNCTN(YltGYl) 
Y2=F(N-1)*DELY+YL0 
CALL     FUNCTN(Y2,GY2) 
16    N  =  N-1 
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IF(N-l)    102,102,103 

103  I F<( GY2-GY1)*XM0DE)    19,19,20 

19  YL0=Y1+DELY 
Y1=Y2 
GY1=GY2 
Y2=F(N-1)*DELY+YLC 

I F(Y2-B)     104, 104, 105 
105    GY2=1.0E+35*XMODE 
GO    TO    16 

104  CALL     FUNCTN(Y2,GY2) 
GO    TO    16 

20  YHI=Y2-DELY 
Y2=Y1 
GY2=GY1 
IF(N-2)    34,34,35 

34  FNM2=0.0 
GO   TO    36 

35  FNM2=F(N-2) 

36  Y1=FNM2*DELY+YL0 
CALL     FUNCTN(Y1,GY1) 
GO    TO    16 

102    IFU  GY2-GY1)*XM0DE)     110,111,111 

110  YYBEST=Y2 
BEST=GY2 
RETURN 

111  YYBEST=Y1 
BEST=GY1 
RETURN 

41  Y1=A 

CALL    FUNCTN<Y1,GY1) 

Y2=B 

CALL    FUNCTN(Y2,GY2) 

GO   TO    102 

42  Y1=A 

CALL     FUNCTN(Y1,GY1) 

GO   TO    111 

END 
CFUNCTN 

SUBROUTINE    FJNCTN(Y,3Y) 

DIMENSION    FM( 10001) ,FNM1 ( 10001 ), 0 NOFX N(l 0001 ) 

REAL     *    8    TAB( 4) 

COMMON    /CON/TAB 

COMMON    XN,N,0N,  YN ,RN, S LQ  W  ,  XHI GH,  D  E  LX, DL0W,DHIGH,  DELD,NUM3ER 

COMMON    YL0H,YHlGH,0ELY,NUFF,FN,JT0P,NNFLAG,NMAX,XLAM 

COMMON    A,B,C,D,E, SUMELT,KQDE 

EQUIVALENCE    ( FN ( 1 ) ,FNM1{ 1  J  ,DN0FXN(1U 

DN=Y 

N0UTPE=6 
C  JUMP    MAY    BE    SET    TO    1     IN    TRANFM    TO     INDICATE    AN    INFEASIBLE    XN-1 

C  JUMP=0 

CALL     TRANFM 

CALL    STGRET 
C  XM0DE=1    FOR    MIN,     -1    FOR    MAX 

C  IF(JUMP.EQ.l) RN=999999*XM0DE 

IF(N-l)     146,146,401 
401    K=(YN-YLOW)/OELY+1.00  1 

IF(K)     212,212,209 
212    WRITE(N0UTPE,210)    N,XN,DN,YN 

210    F0RMAT(///9H    AT    STAGE  15, 9H    HITH   XN=E15.8,8H    AND    DN=E15«8,12H    XN-i 
1EQUALSE16.8,21H    AND    IS    OUT    OF    LIMITS///) 

2002  FORMAT!1!*,'     ERROR    HALT') 

2003  WRITE(N0UTPE,2002) 
STOP 

209    IF(K-JTOP)    205,206,212 

206  RNM1=FNM1(JT3P) 
GO   TO    207 

205    X=K-1 

X=YLOW+X*DELV 

RNM1  =  FNM1<K)«-(FNM1(K«-1I-FNM1(K))*  (YN-XJ/DELY 

207  GO    TO    (21,22,23,241  ,NNFLAG 

21  QN=RN+RNM1 
GO   TO    137 
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22  QN=RN*RNM1 
GO   TO    137 

23  IF(RN-RNMl) 

141  ON=RN 

GO   TO    137 

142  QN=RNM1 
GO   TO    137 

24  IF(RN-RNMl) 

145  QN=RNM1 
GO   TO    137 

146  QN=RN 
137   GY=QN 

RETURN 
END 


141  ,141,142 


145,145,146 
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CSTORE 

SUBROUTINE    STORE 

DIMENSION    FN  (  10QO1)  tFNMK  10001  )f  3  NOFX  N  ( 10001 )  ,  At  1  01)  ,  8(10  1 )  ,C  (  101 ) 

DIMENSION    0(101  )  ,E(101) 

REAL     *    8    TAB(  4) 

COMMON    /C0N/TA8 

COMMON    XN,N,QN,YN,RN,  SLOW ,  XHIGH,  OELX, DL0W»DHI6H» 0  ELD t NUMBER 

COMMON    YLOW,YHIGH ,DEL Y ,NU FF, FN, J T 0 P,MNFL AG VNN AX, X  LAM 

COMMON    AtBfCtCEf  SUMELTf  KCDElf  K30E2 

EQUIVALENCE    (  FN(  1  ) ,  FN  MM  1  )  ,DNOFXM(  D) 

INDATA=8 

N0UTPE=6 

READ     (INDATA,80S)    K0DE1,K0D£2 
808    FORMAT    (213) 

IF(KODEl.EQ.O)     WR  ITE<  NOUTPE  ,  110) 

I  FU0DE1.EQ.1  )    WRITE(NOUTPE,101) 

IF(K0DE1.EQ.2  )  WRITE( NOUTPE ,202) 

IF(K0DE1.EQ.3)     WR  ITE(  NOUTPE  ,  303) 
110    FORMAT    (30H    THE    OBJECTIVE     IS    COST  ) 

101    FORMAT    (30H    THE    OBJECTIVE     IS    SMA  ) 

202    FORMAT     (30H    THE    OBJECTIVE     IS    MS*T  ) 

303    FORMAT    (30H    THE    OBJECTIVE     IS    AVAILABILITY) 

IF(KODE2.EQ.O)     WRITE<  NOUTPE, 900) 

I  FU0DE2.EQ.1  )    WR  IT  E(  NOUTPE  i  901) 

IF<KQDE2.EQ.2 )     WR ITE( NOUT PE ,902) 

IF(K0DE2.  EQ.3  )    WRITE(  NOUTPE,  903) 

900  FORMAT    (30H    CONSTRAINT    ON    COST  ) 

901  FORMAT    (30H    CONSTRAINT    ON    SMA  ) 

902  FORMAT     (30H    CONSTRAINT    ON    MSRT  ) 

903  FORMAT     (30H    CONSTRAINT    ON    AVAILA3ILITY  ) 
DO    333    K=1,NMAX 

READ     (INDATA,999)     A (K ) ,C ( K  ) , E( K)  , B (K)  , D(  K ) 
999    FORMAT    (5F10.4) 

WRITE(NOUTPE,  999)     A  (  K  )  ,B  (  K)  ,  C  (  K)  ,  D  (  K)  ,E(K) 
333    CONTINUE 

SUMELT=0. 

DO    30     I=1,NMAX 
30    SUMELT  =  SUMELT*E  (  I  )*A(  I  )*6<  I) 

WRITE(6,765)SUMELT 
765    FORM ATCF10.3) 

RETURN 

END 
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,STGRET 

SUBROUTINE    STGRET 

DIMENSION    FM  {  10001) ,FNMll 10001 ),} NOFX N( 10001 )  ,  A(  1  01 )  ,  B  (  10  1  ) ,  C  (  101  ) 
DIMENSIGN    0(101  )  ,E<101) 
REAL     *    8    TAB(4) 
COMMON    /CON/TAB 

COMMON    XN,N,0N,YN,RN, SLO W ,XHIGH, 0 E LX, DLQW, DHI GH, D ELD , NUMB ER 
COMMON    YLOWtYHIGH  ,DEL  Y  ,NU  FF,  FN  ,  J  T  OP,  NJNFL  AG,NM  AX,  XLAM 
COMMON    A,B,C,D,E, SUME LT, KCDcl , KOD E2 
EQUIVALENCE    (  FN  (  1  ) ,  FNMH  1 )  ,DNOFXM(  1)) 
I  FU0DE1.EQ.3  )G0    TO    77 
IDN=INT(DN) 
A6=A(N)*B(N) 
I  S=IDN 
TERM=0. 
TEMP=0, 

I F( I S.LT.O)     30    Tj    11 
TERM=EXP(-AB) 
TEMP=TERM 

I FC I S.EQ.O)    30    TO    11 
DO    10     1=1,  IS 
TEMP=TEMP*AB/I 
IF(TERM.GE.. 99999)    GD     TO     11 
C  IF(TEMP. LE-. 00001  )    GO    TO     II 

10  TERM=TERM    +  TEMP 

11  CDF=TERM 
20    CONTINUE 

IFUODEl.NE.l  )G0     TO    40 
C    FOLLOWING    CHANGED    120  133 
C  SMA=  (AB*(  l.-TEMP)  -K  IDN-AB)*(1.-:0F  )  )/ A8 

SMA=(1.-TEMP)*<  IDN-A8)*(  l.-COFJ/AB 

RN=E (N)*AB*S^A/SUMELT 

RETURN 
40    TWUS  =  ( l.-CDF)*< AB*AB-2.*AB*I0N+I3N*(IDN+1)  )/(2.*A(N)  ) 
X*TEMP*B(N)*(AB-IDN)/2. 

IF(K0DE1.EQ.3 )G0     TO    5 0 

RN=E(N)*TWUS/SUMELT 

RETURN 
50   CONTINUE 

AMSRT  =  TWUS/A8 

AMTBF=1./A(N) 

AMTTR=D(N) 

RN=( AMTBF)/(AMTBF+AMTTR+AMSRT) 

RETURN 
77    RN=C (N»*DN 
C  IF(K0DE2.EQ.l .AND.N.EQ.l. AND.YN. LT.XN )RN=999999 

RETURN 

END 
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CDLIMIT 

SUBROUTINE    DLIMIT 

DIMENSION    FN(10OOl),FNMU  10301 ),  D  NOFXNU OOOl )  ,  A(  1  01 J  ,  3(  10  1  ) ,  C  (  10] 

DIMENSION    D(101) ,  E( 10 II 

REAL    *    8    TAB(4) 

COMMON    /CON/TAB 

COMMON    XN,N,3N, YN,RN, SL0W,XHlGH,3ELX,DL0W,DrtIGH,DELD,NUMBER 

COMMON    YLOW, YHI GH, DEL Y,N UFF, FN, J T OP, ^ NFL AG,NM AX, X  LAM 

COMMON    A,B,C,D,E, SUMELT, KCDcl,K3DE2 

EQUIVALENCE    (  FN  (  i  J ,  FNMH  1 )  ,DNOFXN{  1)) 

IF(K0DE2.NE.3 )G0    TO    3 7 

DLOW=0, 

DELD=1. 

DHIGH=AMIN1(XN/C(  N)  ,50.) 

RETURN 
37    AB=A(N)*B(N) 

DO   555    ID=1,30 

IS=ID-1 

TERM=EXP(-A8) 

TEMP=TERM 

I F( I S.EQ.O)    GO    TO    11 

DO    10    1=1  ,IS 

TEMP=TEMP*A3/I 

IF(TERM. GE.. 99999  )    G3    TO    11 
C  IF(TEMP.LE.. 00001  J    GO    TO    11 

10  TERM=TERM    +TEMP 

11  CCF=TERM 

40    TWUS  =  (  U-CDF)  *(  AS  *AB-2.* AB*I  S+  IS*  (  IS*- 1  )  I  /  (  2.  *  A<  N)  J 
X*TEMP*B(N)*(A8-IS  )/2. 

IF(KQDE2.EQ.2.AND.E(N)*TWUS/SUMELT.LE .XN)GO   TO   556 
C      CHANGED    FOLLOWING    TO    555    FROM    556 

IF(K0DE2.EQ.2)G0    TO    555 

IF<K0DE2.EQ.3)G0    TO    554 

SSS= (l.-TEMP)+< IS-AB) *<1.-CJF)M3 

RR=E(N)*AB*SSS/SU*ELT 

Y=XN-RR 

DLOW=0. 
C    THE    NUMBER    IN    THE    NEXT    LINE    HAS    A    BIG    EFFECT.    IT    SHOJLD    NOT. 

IF(N.EQ.1IDL3W=1000. 

IF(N.EC.l.AND.Y.GT.O.  . AN D. ID. LT. 3 0  ) GO    TO    555 

IF(N.EC.1.ANI3.Y.LE.  0.  )DL0W=I5 
136    DHIGH=30.+DL3W 

IFOLQW.EQ.1000.  )  DHIGH=1000. 

DELD=1. 

RETURN 

554  AMSRT=TWUS/A6 
AMTBF=1./A(N) 
AMTTR=D(N) 

AV=AMTBF/(AMTBF+AMTTR+AMSRT) 
IF(AV.GE.XN)     GO    TO    556 

555  CONTINUE 
IS=30, 

556  DLOW=IS 
DELD=1. 
DHIGH=50. 
RETURN 
END 

8L0CK    DATA 
COMMON    /CON/TAB 

REAL     *    8    TA8(4I     /•  SUM1,1       MULT  •  ,  •  M  AXMI  N»  ,  '  MI  NMAX*  / 

END 
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TRANFM 

SUBROUTINE    TRANFM 

DIMENSION    FN  U0001),FNM1(  10001  ),0NOFXN<  10001)  ,  A(  1  01 )  ,  8  ( 10  1 ) ,  C  (  101  ) 

REAL     *    8    TAB(4) 
COMMON    /CON/TAB 

COMMON    XN,N,0N,YN,RN,  SLQ W , XHIGH, D E LX, DLOW , DH I GH, 0  ELD t NUMBER 
COMMON    YL OW, V  HI GH, DEL Y, NU FF, FN, J T OP, M NFL  AG, ^M AX, X  LAM 
COMMON    A,B,C,D,E, SUME LT, KCDE1 , K3D E2 
EQUIVALENCE    (  FN  (  1  ) ,  FNMH  1  )  ,DNOFXM<  1)) 
I F(KODE2.NE.O )G0    TO    12 
YN*XN-C(N)*DN 
RETURN 
12    IDN=INT(DN) 
AB=A(N)*B(N) 
IS=IDN 
TERM=0. 
TEMP  =  0. 

I  F<  I  S.LT.O)    GO    TO    11 
TEPM=EXP(-AB) 
TEMP=TERM 

I F( I S.EQ.O)    GO    TO    11 
DO    10     1=1, IS 
TEMP=TEMP*A8/I 
IF(TERM. GE.. 99999 )    GO    TO    11 
C  IFtTEMP. LE.. 00001  )    GO    TO     11 

10  TERM=TERM    +TEMP 

11  CDF=TERM 

40    TWUS  =  (1.-C0F)  *(  AB*AB-2,*AB*IDN*IDN*(IDN"H)  )/(2.*A(N) ) 
X+TEMP*B(N)*(AB-IDN)/2. 
IF(K0DE2.EQ.2)    GO    TO    50 
IF(K0DE2.EQ.3 )    GO    TO    51 
SMA=(  l.-TEMPJ-M  IDN-A8  )*(  1. -CDF)/ 4  8 
YN=XN-E(N)*AB*SMA/SUMELT 

IF(YN.LT.O)    YN=0.0 

C  IF(N.NE.l.AND.YN.LT.0.OR.N.,NE.l.*ND.YN.GT.l.)  JUMP  =  1 

C  IF(N.E0.1.AND.YN. GT.O )JUMP=1 

C  IF(JUMP.EQ.1)YN=0. 

RETURN 
51    AMSRT=TWUS/AB 

AMTBF=1./A(N) 

AMTTR=D(N) 

AV=( AMTBF)/(AMTBF+AMTTR*AMSRT) 

YN=XN/AV 
C  AV    IS    AVAIL 

RETURN 
50    YN=XN-E(N)*TWUS/SLMELT 

IF(YN.LE.0.)YN=0. 

RETURN 

END 
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DIMENSION    A(  20)  ,8(20)  ,  C(  2  0  ) ,  D(  23  )  ,  E(  2  0  ) ,  MINQ  (  20) 

DIMENSION    T(10,30)  ,F(  1200)  ,DX(13,  1200  )  ,TN<10,31) 

INTEGER    BUDGET, B6 ,BBP 1,DX 

REAL*8    TIME,CTIME 

INDATA  =  8 

N0UTPE=6 

READ( INDATA,888)     NMAX , BUDGET 
888    F0RMAT(2I5> 

READ    (INDATA,777)     (MINQ( I) t 1-1 tMMAXJ 
777    F0RMAT(20I5) 

DO    333    K=1,NMAX 

READ     (INDATA,999)     A (K )  ,C ( K ) , E( K)  ,  B (K)  , D< K ) 
999    FORMAT     (5F10.4) 
C  WRITE(NCUTPE,999)     A (K ) ,B ( K ) , C( K) , D (K )  , E(  K ) 

333    CONTINUE 

TIME=CTIME(1> 

SUMELT=0. 

DO    3  0    1=1, NMAX 
30    SUMELT=SUMELT+E (I )*A< I )*3( I) 


DO    10    1  =  1, NMAX 

DXU  ,1  1  =  0 

DO    110    KP1=1, 31 

K=KP1-1 

AB=A< I  )*B(I) 

TERM=EXP(-A6) 

TEMP=TERM 

IF(K.EQ.O)    G3    TO     11 

DO    14    JI=1,K 

TEMP=TEMP*AB/JI 

c 

IF(TERM.GE.. 999999)    30    TO     11 

c 

IF(TEMP.LE.. 0000001)     GO    TO    11 

14 

TERM=TERM    +TEMP 

11 

COF=TERM 

40    TWUS=(1.-CDF) *( AB*AB-2.*AB*K+K*<  <  +  1 )  )  /  (2.  *AU  )  ) 
X+TEMP*B(I )*(AB-K) /2. 
110    TN(I ,KP1)=E( I )*T*LS/SUMELT 
DO    70    K=l,30 
70    T( I, K)=TN(I,K+1 )-TNU ,K) 
10   CONTINUE 

C  DO   456    IJK=1.10 

C    456    WRITE( 6,555) I TN( K I , I  J K) , KI=1 ,NM4X ) 

C    555    FORMAT(5F12.7) 


DO    80    1=1, NMAX 

LIM=30/MINQ(I ) 

DO    85    K=l ,LIM 

TMOD=0. 

MM=MINC(I ) 

DO    90    J=1,MM 
90    TMOD=TMOD+TlI , ( K- 1 ) *M INQ ( I )+J ) 

T( I, K)=TMOD 
85    CONT  INUE 

C( I) =MINQ(I )*C(  I  ) 
80   CONTINUE 

IBD=BUDGET+1 
DO    100    B8P1=1,IBD 
6B=BBP1-1 
XMIN=0. 
ISTAR=0 
TRY=0. 

DO    210    1=1, NMAX 
IF( BB-C(I) .LT.O)    30    TO    210 
IF(BB-C<n  .GT.O)    GO   TO    211 
TRY=AMIN1( TRY,T(I, 1)1 
GO    TO    212 
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211  TRY=T(I,DX(I,8BP1-C(I  ))  +  l  )+F(BBPl-C(I  )  ) 

212  IF    (TRY.GT.X^IIN)     GO    TO    210 
XMIN=TRY 

I  STAR  =  I 
210    CONT  INUE 

F (BBP1 )=XMIN 

19    DO    21     I»1,NMAX 

IF    (  ISTAR.EQ.O.QR.BB-C(ISTAR)  .LE.O)    DX(I,BBP1)  =  0 
IF    (  ISTAR.EQ.O.QR.BB-C(I  STARJ.LE.  0)    SO   TO    21 
OXU  , BBP1)=DX( 1,3 ePl-Cil STAR)  ) 
21    CONTINUE 

IFUSTAR.EQ.OJGO    TO    100 

IF(BBPl-C(ISTAR).LE.l  )    DX<  ISTAR,  3  B+l)  =1 

IF(BBPl-CUSTAR)  .LE.l)    GO    TO    103 

DX( ISTAR,BBP1 )=AMIN010X< ISTAR,BB? 1-CU  ST  ARI  )  ♦  1 ,  3  0  ) 

100  CONTINUE 
KLCW=IB0-5 
KHIGH=  IBD 

DO    101     KBD=KLOW.KHIGH 

TIME=TIME-CTIME( 1  ) 
C  WRITE(6,3456)TIME 

3456    FORMAT( •TIME=« ,F12.9I 

WRITE(6,1000)(DX( I , KBD ) , I =1 ,NMAX) 

WRITE  (6,1000)  (MINC(I)  ,1  =  1  ,MAX) 
1000    F0RMATU0I5) 

AMSRT=0. 

DO    88     1=1  ,NMAX 
88    AMSRT=AMSRT+TN( I , MINQ ( II *DX( I , KBD ) ♦ U 

KBDM1=KBD-1 

WRITE     (6,2000)     AMSRT,K3DM1 
2000    FORMAT(IX,'MSRT=       SFIO.S,1         BJDGET=       »,I7) 

101  CONTINUE 
RETURN 
END 

FUNCTION    CTIME(  I  ) 
REAL*3    CTIME 
DATA     IFL/O/ 

IF(     IFL.NE.O     )    GO    TO    10 
IFL    =     1 
CTIME    =    O.ODO 
CALL     SETIME 
RETURN 
10    CALL    GETIMEUTIME  ) 

CTIME    =    -ITIME*0«  000026DO 

RETURN 

END 
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